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We propose a mechanism to generate large-scale magnetic fields with correlation lengths of 100 
kpc. Domain walls with QCD-scale internal structure form and coalesce obtaining Hubble-scale cor- 
relations while aligning nucleon spins. Because of strong CP violation, the walls are ferromagnetic, 
which induces electromagnetic fields with Hubble size correlations. The CP violation also induces a 
maximal helicity (Chern-Simons) which supports an inverse cascade, allowing the initial correlations 
to grow to 100 kpc today. We estimate the generated electromagnetic fields in terms of the QCD 
parameters and discuss the effects of the resulting fields. 
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Introduction. The source of cosmic magnetic fields 
with large scale correlations has remained somewhat of 
a mystery There are two possible origins for these 
fields: primordial sources and galactic sources. Primor- 
dial fields are produced in the earlier universe and then 
evolve and are thought to provide seeds which gravita- 
tional dynamos later amplify. Galactic sources would 
produce the fields as well as amplify them. Many mech- 
anisms have been proposed 1^-^, however, most fail to 
convincingly generate fields with large enough correla- 
tion lengths to match the observed microgauss fields with 
~ 100 kpc correlations. We present here a mechanism 
which, although probably requiring a dynamo to produce 
microgauss fields, generates fields with hundred kilopar- 
sec correlations. We present this mechanism as an appli- 
cation of our recent understanding of QCD domain walls, 
which will be described in detail elsewhere Q . 

1. Sometime near the QCD phase transition, Tqcd ~ 
1 GeV, QCD domain walls form. 

2. These domain walls rapidly coalesce until there re- 
mains, on average, one domain wall per Hubble vol- 
ume with Hubble-scale correlations. 

3. Baryons interact with the domain walls and align 
their spins along the domain walls. 

4. The magnetic and electric dipole moments of the 
baryons induce helical magnetic fields correlated 
with the domain wall. 

5. The domain walls decay, leaving a magnetic field. 

6. As the universe expands, an "inverse cascade" 
mechanism transfers energy from small to large 
scale modes, effectively increasing the resulting cor- 
relation length of the observed large scale fields. 

We shall start by discussing the "inverse cascade" mecha- 
nism which seems to be the most efficient mechanism for 
increasing the correlation length of magnetic turbulence. 
After presenting some estimates to show that this mecha- 
nism can indeed generate fields of the observed scales, we 
shall discuss the domain wall mechanism for generating 
the initial fields. 



Evolution of Magnetic Fields. As suggested by Corn- 
wall 1^, discussed by Son [Q and confirmed by Field 
and Carroll ||^ , energy in magnetic fields can undergo an 
apparent "inverse cascade" and be transfered from high 
frequency modes to low frequency modes, thus increas- 
ing the overall correlation length of the field faster than 
the naive scaling by the universe's scale parameter R{T). 
There are two important conditions: turbulence must be 
supported as indicated by a large Reynolds number Re, 
and magnetic helicity (Abelian Chcrn-Simons number) 
H — / A • Bd'^a; is approximately conserved. The im- 
portance of helicity was originally demonstrated by Pou- 
quet and collaborators |^]. The mechanism is thus: the 
small scale modes dissipate, but the conservation of helic- 
ity requires that the helicity be transfered to larger scale 
modes. Some energy is transfered along with the helicity 
and hence energy is transported from the small to large 
scale modes. This is the inverse cascade. The reader is 
referred to |^-^ for a more complete discussion. 

In the early universe. Re is very large and supports tur- 
bulence. This drops to Re w 1 at the e^e~ annihilation 
epoch. To « 100 eV After this point (and throughout 
the matter dominated phase) we assume that the fields 
are "frozen in" and that the correlation length expands 
as R while the field strength decays as R^^ . Note that 
the inverse cascade is only supported during the radiation 
dominated phase of the universe. 

Under the assumption that the field is maximally he- 
lical, these conditions imply the following relationships 
between the initial field i?rms(7i) with initial correlation 
l{Ti) and present fields today (Tnow ~ 2 x 10^* eV) 
-BrmslTnow) with Correlation l{T^ow) 
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As pointed out in , the only way to generate turbu- 
lence is either by a phase transition Ti or by gravitational 
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instabilities. We consider the former source. As we shall 
show, our mechanism generates Hubble size correlations 
li at a phase transition T^. In the radiation dominated 
epoch, the Hubble size scales as T~^. Combining this 

with (H), we see that Znow oc T~^^^; thus, the earher the 
phase transition, the smaller the possible correlations. 

The last phase transition is the QCD transition, Ti = 
Tqcd « 0.2 GeV with Hubble size IITqcd) ~ 30 km. 
We calculate (^) the initial magnetic field strength to be 
B,^s{T,) « eA|c^/(eAQCD) ~ (10i^G)/(eAQCD) where 
^ is a correlation length that depends on the dynamics 
of the system as discussed below and Aqcd ~ 0.2 GeV. 
With these estimates, we see that 



it couples to the anomaly and is the dominant player in 
aligning the magnetic fields. The potential 
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today. One could consider the electroweak transition 
which might produce 100 pc correlations today, but 
this presupposes a mechanism for generating fields with 
Hubble-scale correlations. Such a mechanism does not 
appear to be possible in the Standard Model. Instead, 
the fields produced are correlated at the scale which 
can produce only ^ 1 km correlations today. 

These are crude estimates, and galactic dynamos likely 
amplify these fields. The important point is that we can 
generate easily the 100 kpc correlations observed today 
provided that the fields were initially of Hubble size cor- 
relation. Unless another mechanism for amplifying the 
correlations of magnetic fields is discovered, we suggest 
that, in order to obtain microgauss fields with 100 kpc 
correlation lengths, helical fields must be generated with 
Hubble-scale correlations near or slightly after the QCD 
phase transition Tqcd- The same conclusion regarding 
the relevance of the QCD scale for this problem was also 
reached in . The rest of this work presents a mecha- 
nism that can provide the desired Hubble size fields, jus- 
tifying the estimate (^) . We shall explain the mechanism 
and give simple estimates here. See [|| for details. 

Magnetic field generation mechanism. The key play- 
ers in our mechanism are domain walls formed at the 
QCD phase transition that possess an internal structure 
of QCD scale. We shall present a full exposition of these 
walls in 1^ but to be concrete, we shall discuss an axion 
wall similar to that described by Huang and Sikivie ^ . 

We start with a similar effective Lagrangian to that 
used by Huang and Sikivie except we included the effects 
of the rj' singlet field which they neglected: 



-|' + ^Tr|a^U|2-y(U,a) (4) 



where d = f^^a is the dimensionless axion field and the 
matrix U = exp(i77' -I- iir^ X^) contains the pion and rj' 
fields (to simplify the calculations, we consider only the 
SU{2) flavor group). Although the rj' field is not light. 
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was first introduced in It should be realized that 
iln(det(U)) = iln(det(U)) + 27rn is a multivalued func- 
tion and we must choose the minimum valued branch. 
Details about this potential are discussed in but sev- 
eral points will be made here. All dimensionful parame- 
ters are expressed in terms of the QCD chiral and gluon 
vacuum condensates, and are well known numerically: 
M = -diag(m^|(q'g*)|) and E = {bas/{32Tr)G^). This 
potential correctly reproduces the Veneziano-Witten ef- 
fective chiral Lagrangian in the large Nc limit ||l^ ; it re- 
produces the anomalous conformal and chiral Ward iden- 
tities of QCD; and it reproduces the known dependence 
in 6 for small angles ||l^. We should also remark that 
the qualitative results do not depend on the exact form 
of the potential: domain walls form naturally because of 
the discrete nature of the symmetries [^,|ll|,^ . 

The result is that two different types of axion domain 
walls formJ0|. One is almost identical to the one dis- 
cussed in [Hi with small corrections due to the 77'. We 
shall call this the axion/pion (a^) domain wall. The sec- 
ond type, which we shall call the axion/eta' {ari>) domain 
wall is a new solution characterized by a transition in 
both the axion and rj' fields. The boundary conditions 
(vacuum states) for this wall are a(— 00) = r)'(— 00) = 
and a{oo) = fj'{oo) = ±7r with tt*' = at both bound- 
aries. The main difference between the structures of the 
two walls is that, whereas the domain wall has struc- 
ture only on the huge scale of rn~^, the rj' transition in 
a,f has a scale of m~,^ ~ ■'^qcd- "^^^ reason is that, in the 
presence of the non-zero axion (9) field, the pion becomes 
effectively massless due to its Goldstone nature. The rj' 
is not sensitive to 9 and so its mass never becomes zero. 
It is crucial that the walls have a structure of scale Aq^j-, : 
there is no way for the wall to trap nucleons because 
of the huge difference in scales but the a,,/ wall has ex- 
actly this structure and can therefore efficiently align the 
nucleons. 

The model we propose is this: Immediately after the 
phase transition, the universe is filled with domain walls 
of scale Jqcq . As the temperature drops, these domain 
walls coalesce, resulting in an average of one domain 
wall per Hubble volume with Hubble-scale correlations 
|]ll| , p^ . It is these a^/ domain walls which align the dipole 
moments of the nucleons producing the seed fields. 

The following steps are crucial for this phenomenon: 
1) The coalescing of QCD domain wall gives the fields 
TT^ , ri' Hubble-scale correlations. 2) These fields interact 
with the nucleons producing Hubble-scale correlations of 
nucleon spins residing in the vicinity of the domain wall. 
(The spins align perpendicular to the wall surface.) 3) 
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Finally, the nucleons, which carry electric and magnetic 
moments (due to strong CP violation), induce Hubble- 
scale correlated magnetic and electric fields. 4) These 
magnetic and electric fields eventually induce a nonzero 
helicity which has the same correlation. This helicity 
enables the inverse cascade. 

Quantitative Estimates. As outlined below, we have 
estimated the strengths of the induced fields in terms 
of the QCD parameters We consider two types of 
interactions. First, the nucleons align with the domain 
wall. Here we assume that the fluctuations in the nucleon 
field \[' are rapid and that these effects cancel leaving the 
classical domain wall background unaltered. Thus, we 
are able to estimate many mean values correlated on a 
large scale on the domain walls such as {^j5(Txy'^) and 
(*7275*) through the interaction i'(i^-mAre''''(^''^'s)^'. 

To estimate the magnetization of the domain wall, we 
make the approximation that the wall is flat compared to 
the lengths scales of the nucleon interactions. By assum- 
ing that momentum is conserved in the wall, we reduce 
our problem to an effective 1 + 1 dimensional theory (in 
z and t) which allows us to compute easily various mean 
value using a bosonization trick ||lj,|l^. The result for 
the mean value (^75cra;y\l/) for example is [^: 
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where /i ~ ttiat is a dimensional parameter originating 
from the bosonization procedure of the corresponding 2D 
system and the parameter Aqqq ~ J dkxdky comes from 
counting the nucleon degeneracy in the x-y plane of a 
Fermi gas at temperature Tc ~ Aqcd- These mean values 
are only nonzero within a distance AqJ^j-, of the domain 
wall and are correlated on the same Hubble-scale as the 
domain wall. 

From now on we treat the expectation value ^ as 
a background classical field correlated on the Hubble- 
scale. Once these sources are known, one could calculate 
the generated electromagnetic field by solving Maxwell's 
equations with the interaction 

Ant = ^(d**cr^..75* + M**«crM^*)^P^ + ^(iD)^^' (7) 

where d^j, (/^*) is effective electric (magnetic) dipole mo- 
ments of the field 5*. Due to the CP violation (nonzero 
6) along the axion domain wall, the anomalous nucleon 
dipole moment in (|^) d>p ~ ^\si ~ is also nonzero [ p^ . 
This is an important point: if no anomalous moments 
were induced, then only charged particles could generate 
the magnetic field: the walls would be diamagnetic not 
ferromagnetic as argued in JlTf and Landau levels would 
exactly cancel the field generated by the dipolcs. 

Solving the complete set of Maxwell's equations, how- 
ever, is extremely difficult. Instead, we use simple di- 
mensional arguments. For a small planar region of area 
S,'^ filled with aligned dipoles with constant density, we 



know that the net magnetic field is proportional to 
since the dipole fields tend to cancel, thus for a flat sec- 
tion of our domain wall, the field would be suppressed 
by a factor of (C^qcd)"^- For a perfectly flat, inflnite 
domain wall (^ —^ oo), there would be no net fleld as 
pointed out in ||l^. However, our domain walls are far 
from flat. Indeed, they have many wiggles and high fre- 
quency modes, thus, the size of the flat regions where the 
flelds are suppressed is governed by a correlation ^ which 
describes the curvature of the wall. Thus, the average 
electric and magnetic flelds produced by the domain wall 
are of the order 
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^Aqcd 



[d* (^'crp^75^') + Ai* (^icT^,.*)] (8) 



where ^ is an effective correlation length related to the 
size of the dominant high frequency modes. 

To estimate what effective scale ^ has, however, re- 
quires an understanding of the dynamics of the domain 
walls. Initially, the domain walls are correlated with 
a scale of Aq^.^. As the temperature cools, the walls 
smooth out and the lower bound ^_ (i) for the scale of the 
walls correlations increases from ^_(0) ~ Aqqj-,. This in- 
crease is a dynamical feature, however, and is thus slow. 
In addition, the walls coalesce and become correlated 
on the Hubble-scale generating large scale correlations. 
Thus the wall has correlations from ^_ (t) up to the upper 
limit set by the Hubble-scale. We expect that ^ ^ Hub- 
ble size at the time that the fields are aligned and that 
the suppression is not nearly as great as implied in [ pT[ . 
Note that, even though the effects are confined to the re- 
gion close to the wall, the domain walls are moving and 
twisted so that the effects occur throughout the entire 
Hubble volume. 

The picture is thus that fields of strength 
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(9) 



are generated with short correlations ^, but then domains 
are correlated on a large scale by the Hubble-scale modes 
of the coalescing domain walls. Thus, strong turbulence 
is generated with correlations that run from Aqcd up to 
the Hubble-scale. 

Finally, we note that this turbulence should be highly 
helical. This helicity arises from the fact that both elec- 
tric and magnetic fields are correlated together along the 
entire domain wall, (E) ^ (A)/r where (A) is the vec- 
tor potential and t is a relevant timescale for the electric 
field to be screened (we expect t ^qcd as we discuss 
below). The magnetic helicity density is thus 
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Note carefully what happens here: The total helicity was 
zero in the quark-gluon-plasma phase and remains zero in 
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the whole universe, but the helicity is separated so that in 
one Hubble volume, the helicity has the same sign. The 
reason for this is that, as the domain walls coalesce, ini- 
tial perturbations cause either a soliton or an antisoliton 
to dominate and fill one Hubble volume. In the neighbor- 
ing volume, there will be other solitons and antisolitons 
so that there is an equal number of both, but they are 
spatially separated which prevents them from annihilat- 
ing. This is similar to how a particle and antiparticle may 
be created and then separated so they do not annihilate. 
In any case, the helicity is a pseudoscalar and thus has 
the same sign along the domain wall: The entire Hubble 
volume has helicity of the same sign. This is the origin of 
the Hubble-scale correlations in the helicity and in B^. 
The correlation parameter ^ which affects the magnitude 
of the fields plays no role in disturbing this correlation. 

Eventually, the electric field will be screened. The time 
scale for this is set by the plasma frequency for the elec- 
trons (protons will screen much more slowly) LUp ^ Aqcd ■ 
The nucleons, however, also align on a similar timescale 
^QCD' helicity is generated on this scale too, 

so the electric screening will not qualitatively affect the 
mechanism. Finally, we note that the turbulence requires 
a seed which remains in a local region for a timescale 
set by the conductivity a cT j Aqcd where for 
T = 100 MeV, c w 0.07 (Ts) and is smaller for higher 
T. Thus, even if the domain walls move at the speed of 
light (due to vibrations), there is still time to generate 
turbulence. 

For this mechanism to work and not violate current 
observations, it seems that the domain walls must even- 
tually decay. Several mechanisms have been discussed 
for the decay of axion domain walls and the 

timescales for these decays are much larger than AqJ^^, 
ie. long enough to generate these fields but short enough 
to avoid cosmological problems. QCD domain walls ||] 
are quasistable and may nicely solve this problem. We as- 
sume that some mechanism exists to resolve the domain 
wall problem in an appropriate timescale. Thus, all the 
relevant timescales are of the order Aq^j-, except for the 
lifetime of the walls and thus, although the discussed in- 
teractions will affect the quantitative results, they will 
not affect the mechanism or substantially change the or- 
der of the effects. 

Conclusion. We have shown that this mechanism can 
generate the magnetic fields (^ with large correlations, 
though galactic dynamos should still play an important 
amplification role. It seems that the crucial conditions 
for the dynamo to take place are fields B > 10~^° G 
with large (100 kpc) correlations. From (m we see that 
we have a huge interval 1 < ^Aqcd <C 10 of ^ to seed 
these dynamos. Also, if ^ is small, then this mechanism 
may generate measurable extra-galactic fields. 

We mention two new points that distinguish this mech- 
anism from previous proposals [ pO{ . First, the key nu- 
cleon is the neutron which generates the fields due to an 



anomalous dipole moment induced by the CP violating 
domain walls. The nucleons thus make the wall ferro- 
magnetic, not diamagnetic as discussed in ||l^. Second, 
the interaction between the domain walls and nucleons 
are substantial because of the QCD scale of the 77' tran- 
sition. There is no way that axion domain walls with 
scales iTT-a^ can efficiently align nucleons at a temper- 
ature Tqcd- 

We should also note that the magnitudes of the fields 
generated by this mechanism are small enough to sat- 
isfy the constraints placed by nucleosynthesis and CMB 
distortions. Thus, domain walls at the QCD phase tran- 
sition, in particular those described in |^ , provide a nice 
method of generating magnetic fields on 100 kpc correla- 
tions today (|^). 
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